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1.  Introduction 

Equations  of  the  Schrodinger  type  arise  in  many  disciplines,  such  as  quantum  mechanics,  fluid 
mechanics,  plasma  physics,  laser  propagation,  acoustics  and  optics  [4,  8, 19,  23,  28,  27,  34,  38,  36, 

37] .  This  paper  is  primarily  concerned  with  the  numerical  solution  of  the  equation 

u, «  riu*,,  (1.1) 

with  Am  a  +  it  and  a  >  0,  and  its  extension  to  higher  dimensions: 

m 

(1.2) 

Jatl 

where  Ai  *  oj  +  *8j,  with  ai  >  0  and  8j  real,  and  af  +  Ilf  yf  0  for  /  ■  !,••*, m. 

Equation  (1.1)  includes  both  the  heat  equation  ««  ■  iun  and  the  Schrodinger  equation 
«t  v  .  It  is  well  known  that  a  rather  complete  collection  of  stability  results  for  difference 
schemes  exist  for  the  heat  equation  [29]: 

«iaauH.  (1.3) 

and  for  the  advection-di {fusion  equation  [6,  9]: 

«(■«««  +  tu,.  (1.4) 

It  is  our  intention  to  provide  a  similar  collection  of  results  for  the  Schrodinger  equation.  We  propose 
a  collection  of  finite  difference  schemes,  and  analyse  their  accuracy  and  stability  properties.  Some 
of  the  schemes  are  analogous  to  well-known  schemes  for  the  wave  equation  and  the  heat  equation 
but  others  are  unique  to  the  Schrodinger  equation.  This  includes  explicit  and  implicit  schemes, 
2-level  and  3-level  schemes  and  real  and  complex  schemes.  Many  of  these  are  analogous  to  classical 
schemes  for  the  heat  equation  but  some  schemes  are  unique  to  the  Schrodinger  equation.  Von 
Neumann  type  stability  results  are  derived  for  all  the  schemes  and  extensions  to  higher  dimensions 
are  derived  in  most  cases.  Many  of  stability  results  are  quite  different  from  the  corresponding 
results  for  the  heat  equation  and  the  wave  equation. 

The  existence,  uniqueness  and  regularity  properties  of  equations  of  the  Schrodinger  type  have 
been  investigated  in  recent  years  [2,  3,  14,  17,  28,  34,  39].  We  are  mainly  going  to  discuss  finite 
difference  methods  for  such  equations.  Among  the  numerical  methods  for  these  equations,  the 
finite  difference  method  is  not  only  a  basic  one,  but  also  one  of  the  most  extensively  used.  Since 
many  conventional  explicit  schemes  are  unconditionally  unstable  for  the  Schrodinger’s  equation 
[1,  16,  19,  23],  implicit  schemes  have  been  the  most  popular  —  especially  the  Crank-Nieolson 
scheme.  These  results  can  be  found  in  (1,  10,  II,  18,  16,  19,  23,  27].  Recently,  it  has  been  found 
that  stable  explicit  schemes  for  equations  of  the  Schrodinger  type  can  be  derived  if  appropriate 
dissipative  terms  are  added  [8]  and  some  of  these  explicit  schemes  have  been  applied  to  underwater 
acoustics  problems  [7].  The  articles  [16,  40]  investigate  the  existence  and  convergence  of  implicit 
schemes.  In  recent  years,  the  trend  of  applying  spectral  methods  [4,  8,  12,  13,  18,  28,  38,  36,  37, 

38]  and  the  finite  element  methods  [11, 16,  31]  is  increasing  as  well.  Methods  of  lines  methods  have 
also  been  used  [22,  24].  Since  the  solutions  of  nonlinear  equations  of  the  Schrodinger  type  often 
possess  conservation  laws,  attempts  have  also  made  to  construct  schemes  which  satisfy  discrete 
conservation  laws.  To  achieve  this,  M.  Delfour  et.  ai  [11]  modified  the  Crank- Nicolson  scheme 
and  J.M.  Sanz-Serna,  and  V.S.  Manoranjan  [30]  used  the  Leap-frog  technique.  Among  three  level 


schemes,  the  Leap-frog  scheme  hss  been  suggested  [31].  In  [15],  predictor-corrector  schemes  are 
discussed.  For  applications  of  some  of  the  schemes  proposed  b  this  paper,  the  reader  is  referred 
to  the  references  (7,  20,  21,  23]. 

In  Sec.  2,  our  definition  of  stability  b  given,  as  well  as  a  general  method  for  deriving  stability 
results  from  the  characteristic  polynomial  of  a  numerical  scheme.  In  Sec.  3,  the  special  case  o  =  0 
b  considered  first  and  the  issue  of  the  existence  of  stable  2-level  explicit  schemes  b  addressed. 
In  Sec.  4,  we  consider  a  general  two-level  scheme  for  (1.1).  In  Sec.  5,  6  and  7,  we  consider 
3-level  schemes:  the  Leap-frog  scheme,  the  Do-Fort  Frankel  scheme  and  the  backward  difference 
scheme.  In  Sec.  8,  we  separate  the  real  and  imagbary  parts  of  (1.1),  and  consider  schemes  that 
are  specifically  designed  for  the  resulting  system  of  real  equations. 

Throughout  this  paper,  k  denote  the  temporal  mesh  size,  h  the  spatial  mesh  size,  r£  p,  ij  * 
4sb*  |,  0  <  #  <  2r  and  *y  *  nj. 

2.  Definition  of  Stability  and  The  Schnr-Cohn  Theory 

The  usual  von  Neumann  type  definition  of  stability  requires  that  the  roots  Rj  of  the  charac¬ 
teristic  polynomial  of  a  numerical  scheme  satisfy 

\Ri\Zl  +  0(k)  (2.1) 

[29].  While  thb  b  the  appropriate  definition  of  stability  for  provbg  convergence  as  k  and  h  tend  to 
zero,  b  conjunction  with  the  Lax-Equivalence  theorem  [29],  for  practical  computations  with  fixed 
k  and  k  thb  definition  allows  the  numerical  solution  to  grow  with  the  number  of  time  steps  taken. 
For  equations  the  solutions  of  which  are  known  to  be  nonmcreasing  b  time,  as  b  the  case  for  (1.1), 
thb  b  often  undesirable.  Hence,  for  the  stability  of  a  numerical  scheme,  we  shall  require  that  the 
numerical  solutions  also  do  not  grow  b  time.  Tins  definition  of  stability  b  sometimes  known  as 
the  practical  liability  criteria  [29]  and  b  slightly  stricter  than  the  definition  (2.1).  In  the  absence 
of  lower  order  terms  (e.g.  u,  or  u)  or  b  the  limit  as  k  and  h  tend  to  zero,  the  difference  between 
the  two  definitions  of  stability  b  usually  very  slight.  In  what  follows,  we  shall  make  this  definition 
of  stability  more  precise  and  outline  a  procedure  for  systematically  deriving  stability  conditions 
for  a  given  numerical  scheme. 

We  shall  follow  the  methodology  developed  b  [6, 9].  We  define  two  classes  of  polynomials: 
Definition  2.1.  : 

We  shall  call  polynomials  d(g)  with  roots  Rj  Schur  Polynomial*  if  |/i/|  <  1  Vj,  and  Simple 
von  Neumann  Polynomial*  if  |A;-|  <  1  V;  and  the  roots  with  magnitude  equal  to  one  are  dbtinct. 

Let  d(z)  be  the  characteristic  polynomial  corresponding  to  a  particular  scheme,  obtamed  via 
Fourier  analysb  [29]. 

Definition  23.  : 

A  numerical  scheme  is  defined  to  be  itable  if  its  characteristic  polynomial  b  Simple  von  Neu¬ 
mann. 

To  determine  whether  a  polynomial  b  a  Simple  von  Neumann  polynomial,  we  shall  use  the 
theory  of  Schur  [26, 33,  32].  Thb  theory  enables  one  to  determine  conditions  on  the  coefficients  of 
the  characteristic  polynomial  for  it  to  be  Simple  von  Neumann. 

Given  a  polynomial 

¥ 

d(s)  -  flo  +  «l*  +  —  + 


of  degree  v  (with  a„  #  0,  a0  /  0),  one  can  associate  with  4  another  polynomial  4*,  defined  by 


+*(*)  * 

i*=® 

where  a  denotes  the  complex  conjugate  of  a.  The  reduced  polynomial  4i  is  defined  bjr 

By  defintion,  the  degree  of  4i  is  one  less  than  that  of  4-  The  main  results  that  we  need  are  contained 
in  the  following  two  theorems: 

Theorem  14.  4  **  »  Schur  Polynomial  iff  |d*(0)|  >  |d(0)|  and  4i  i*  a  Schur  Polynomial. 

Theorem  1.1.  4  »  •  Simple  von  Neumann  Polynomial  iff  either  (1)  |d*(0)|  >  |d(0)|  and  4i  is  a 
Simple  von  Neumann  Polynomial,  or  (2)  4i  s  0  and  4'  **  •  Schur  Polynomial  (4'  denotes  the 
derivative  of  4  with  respect  to  its  dependent  variable). 

By  repeated  applications  of  the  above  two  theorems,  it  is  possible  to  reduce  the  question  of 
whether  an  n*th  degree  polynomial  is  a  Simple  von  Neumann  Polynomial  to  that  for  a  first  degree 
polynomial,  which  can  be  solved  more  easily  by  analytical  means.  This  procedure  turns  out  to 
be  very  effective  for  determining  stability  limits  of  general  numerical  schemes,  as  compared  to 
first  finding  the  roots  of  the  characteristic  polynomial  explicitly  and  then  determining  conditions 
for  their  absolute  values  to  be  less  than  unity.  Furthermore,  this  last  approach  may  not  even  be 
applicable  for  polynomials  of  higher  degrees  which  arise  in  the  analysis  of  multi-level  schemes  and 
systems  of  equations  (see  Sec.  8).  Finally,  it  is  worth  noting  that  this  reduction  process  preserves 
the  necessity  and  sufficiency  of  the  stability  conditions. 

9.  Stable  Explicit  Schema  for  ut  «=  iuM 

We  shall  first  consider  the  more  special  equation 

til  *  iu„.  (3.1) 

Even  in  this  simple  case,  the  stability  properties  of  some  popular  schemes  are  quite  different  from 
that  for  the  superficially  similar  heat  equation  u<  =  u„. 

S4.  Taylor  Stria  Schema 

Consider  the  basic  Euler  Scheme: 

—  U* 

— — £ — -  *=  iDjV  ,  (3.2) 

u*,.  —  2u*  +  u*  . 

where  Dfv  s  ^ - — .  (3.3) 

The  truncation  error  is  0(4,  h2).  The  corresponding  characteristic  polynomial  is 

d(*)  -  *  -  (1  -  #7).  (3.4) 

Since  the  only  root  is  R  *  1  -if  and  |/2|2  =  1  +  7*  >  1  for  7  /  0,  4(z)  is  not  Simple  von  Neumann, 
and  thus  the  Euler  scheme  is  unconditionally  unstable,  as  is  well-known  |8,  23). 


Hie  Euler  scheme  is  the  first  member  of  s  class  of  schemes  derived  from  the  following  Taylor 
series  expansion  of  a  solution  u(x,  t)  of  (1.1): 

fc*  k* 

u(x,t  +  k)  m  u(x,t)  +  kut  +  -j-a«  +  +  ••■ 

»  u(x,  t )  +  k(iuMM)  +  +  •  •  •  (3.5) 

A  family  of  schemes  can  be  derived  from  (3.6)  by  replacing  the  spatial  derivatives  by  difference 
operators,  namely: 

«•+*  -  «•  '  tl-i  . 

lm  1 

The  truncation  error  is  0(k*,ks).  The  Eluler  scheme  corresponds  to  p  ■  1.  By  retaining  the  next 
term  in  the  expansion  (p  «  2),  we  obtain  the  following  Lax-Wendroff  type  scheme: 

j  (!$)»«. 

The  truncation  error  is  0(4*,**).  The  corresponding  characteristic  polynomial  is 

*(*)■*“  (l  -  *7  -  y)  - 

Since  the  only  root  is  R  •  (l  -  -  »7  and 

\R\7  *  (l  -  y)  +TJ“1  +  7>1  for  TT^O, 

this  scheme  is  also  unconditionally  unstable  [8]. 

Since  the  first  two  members  of  this  family  of  schemes  are  unconditionally  unstable,  it  b 
of  theoretical  interest  to  determine  whether  there  b  any  member  that  b  conditionally  stable. 
Moreover,  thb  has  practical  implications  as  well,  because  simple  explicit  schemes  are  easier  to 
implement  (and  vectorise)  than  implicit  schemes,  especially  for  higher  dimensional  problems. 

The  third  order  (in  time)  scheme  derived  thb  way  has  a  characteristic  polynomial  given  by 


the  condition  that  d(z)  be  Simple  von  Neumann  reduces  to 


which  is  satisfied  if  and  only  if  i  <  v/3.  The  stability  condition  is  thus 


Fot  the  fourth  order  scheme,  we  have 


4>{z) 


7*  *7* 

17  -  -Jr  +  -r  + 
2  6 


After  a  similar  computation,  it  can  be  verified  that  this  scheme  is  stable  if  and  only  if  r  <  . 

We  thus  see  that  there  are  stable  members  in  this  family  of  schemes.  Unfortunately,  while 
the  third  and  fourth  order  schemes  am  conditionally  stable,  they  am  not  very  practical  for  initial 
boundary  value  problems  because  their  spatial  stencils  extends  over  7  (resp.  9)  points,  which  makes 
the  construction  of  stable  numerical  boundary  conditions  difficult. 


This  leads  to  a  natural  question: 

Doe*  there  exist  stable  explicit  scheme*  with  smaller  stencils  9 


U.  Schamas  With  Artificial  Dissipation 

The  answer  to  this  question  is  positive.  As  shown  in  [8],  stable  explicit  schemes  with  a  5-point 
stencil  can  be  derived  from  the  Euler  scheme  by  adding  appropriate  dissipative  terms.  Although 
the  addition  of  dissipation  to  stabilise  a  numerical  scheme  is  rather  natural,  the  question  remains 
as  to  whether  this  is  the  only  way  to  obtain  a  stable  scheme  for  (3.1).  In  what  follows  we  shall 
answer  this  in  the  positive  by  deriving  the  dissipative  schemes  in  [8]  from  a  general  sheme  with  a 
5-point  spatial  stencil  that  satisfies  certain  symmetry  conditions. 

Suppose  an  explicit  scheme  to  solve  (3.1)  takes  the  form 


—  u*t 
'  k  ’ 

1 


(3.6) 


Of  course,  52  cj+tu/+i  must  be  a  consistent  approximation  of  <uM  for  (3.6)  to  be  consistent  with 
(3.1).  Using  Taylor  expansion,  we  find  that  the  following  consistency  conditions  must  be  satisfied: 


**■+* *  0 
1 


■  0 


(3.7) 


If  we  look  for  schemes  with  symmetry,  i.e.  cj+t  «  c;_i,  then  (3.7)  becomes 

W  *°* 

l>0 

2Yll*ei+t  “ 

(X> 


2i 
h * 


(3.8) 
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If  we  take  /  =  -1,0, 1,  i-e.,  using  3  points  to  approximate  tuss,  from  (3.8)  we  find  eo  =  -jt, 
Cj  «  c_i  =  -fr.  This  is  just  the  Euler  scheme  (3.2)  and  is  unconditionally  unstable. 

Now,  let  us  take  /  *»  -2, -1,0, 1,2,  i.e.,  using  5  points  to  approximate  iuM.  Following  the 
same  procedure  as  above,  we  can  easily  find 

2i 

p  • 

where  cz  b  a  constant  to  be  determined.  We  want  to  determine  cz  such  that  scheme  (3.6)  is  stable. 
Writing  cj  «  p,  where  c  *  o  +  iff  with  real  constants  a,  ff,  (3.6)  becomes 

„•+»  •  j 

— L  **  K®  +  W  (U>+J  +  «i-j)  +  (•  “  4(°  +  W)  («*+i  +  «*-»)  («(«  +  »0)  “  2») «"]  • 

(3.9) 

After  some  rearrangements,  scheme  (3.9)  can  be  rewritten  into  another  form: 


(3.10) 


The  last  term  can  be  viewed  as  a  dissipative  term  which  is  added  to  the  unstable  scheme 

<  _  .  fu%i  ~  2“>  +  «i-i  \ 


*  v  ** 


and  the  truncation  error  is  0(k,h9).  The  root  of  the  characteristic  polynomial  for  (3.9)  is 

R  *  1  -  «'ri|  +  (a  +  iff)rrf*. 

Thus  we  have  .  .  .  . 

|ff|*-(l  +  or,*)*  +  (/fn,*-n,)* 

=  1  +  (r*  +  2ra)»j*  -  2^^ffl|•  +  r’(o*  +  ff*)if4. 

The  condition  |f2p  <  1  reduces  to 

(r  +  2o)  -  2 rffi|  +  r(o*  +  02)ni  ^  0, 


from  which  it  follows  that  the  condition  on  r  is 


'■Sffa) 


Or J»fs  +  (fff|  -  l)s  ' 


For  a  finite  stability  interval,  we  must  have  o  <  0.  By  differentiating  it  can  be  verified 
that  g(tf)  cannot  achieve  its  minimum  within  the  interval  0  <  if  <  4.  Thus  the  conditions  on  r 
reduces  to 

r  <  min(f(0),#(4)), 


(3.11) 
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To  obtain  the  scheme  with  the  least  restrictive  stability  condition,  one  should  make  the  right 
hand  side  of  (3.11)  as  large  as  possible.  Obviously,  one  should  take  0  *  1.  Condition  (3.11)  then 
reduces  to 

r  <  min 

The  right  hand  side  achieves  its  maximum  when  -2a  =  -  jh,  or  a  =  — j.  Thus  we  have  for 
the  stability  condition: 

r<\.  (3.12) 

We  summarize  the  results  in  the  following  theorem. 

Theorem  8.1.  For  a ay  real  0  tad  a  <  0,  the  sheae  (3.10)  is  conditionally  stable,  the  necessary  and 
sufficient  condition  being  (3.11).  The  least  restrictive  stability  condition  is  (3.12)  and  is  obtained 
when  a  *  —  \,0  — 


We  can  also  consider  a  similar  dissipative  scheme: 


u*+i  _  u» 

i  i 


i(2 


U"-H  ~  2“"  +  U"-l 


A* 


)  +  (®  +  *P)k  (— 


M"+a  -  <«*4t  j-  6uy  -  4u"-i  + 


)  ,  (3.13) 


whose  dissipative  term  is  different  from  (3.10).  For  (3.13),  the  following  theorem  has  been  proved 
in  (8): 


Theorem  8.8.  The  scheme  (3.13)  is  stable  if  and  only  if  a  <  except  for  the  half  line 
{or  * -\,0  <  0}  ,  and 

^  V'-o  (2a2  +  20s  +  a) 

4(a»  +  ^) 


(3.14) 


The  least  restrictive  stability  constraint  is 


(3.15) 


and  is  obtained  when  a  *  0  «  |. 

U.  Stable  Dissipative  Schemes  In  Mottl-Dlmenslons 
We  next  consider  the  multi-dimensional  equation: 

m 

Ut  *  £ 

1*1 


We  assume  that  the  kj’s  have  the  same  sign. 

Without  loss  of  generality,  we  assume  6;  >  0  (1  *  l,  •  •  ■  ,m).  We  consider  the  natural  exten¬ 
sion  of  scheme  (3.10): 


~ T  L  »£*<  +  (® +  W?  (^»)* «]  •  (316) 

1*1 

Here  j  represents  a  multi-index  (ji,  -  -  -  ,jin),  Dfj  is  the  second  order  centered  difference  operator 
with  respect  to  ji  and  hi  is  the  corresponding  mesh  site. 


—  j 


\ 


Theorem  &A.  Scheme  (3 .16)  is  stable  if  and  only  if  a  <  0  and 


k  <  mis 


2a 
*«  .  * 


IS*  ^+tf-»«S*- 

The  least  restrictive  stability  constraint  is 


k  < 


*£* 


and  is  obtained  when  a  =  and 
Proof.  See  Appendix. 
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4.  A  General  Two-Level  Scheme 


4.1.  The  One  Dimensional  Case: 

We  now  return  to  the  attention  of  the  more  general  equation  (1.1).  We  consider  the  following 
finite  difference  scheme  for  (1.1): 


u"+l  -  u" 


£  («#}  -  + *£.■) + <«■+.  -  *«• + «?-.)  •  <4» 


where  p  is  a  parameter  such  that  0  <  p  <  1.  It  is  easy  to  see  that  when  p «  0,  (4.1)  is  the  explicit 
Euler  scheme,  when  p  =  §,  (4.1)  reduces  to  the  Crank-Nicolson  scheme  and  when  p  *  1,  (4.1)  is 
the  fully  implicit  Backward  Euler  scheme.  The  truncation  error  is  0(4,  h*)  except  it  is  0(4* ,  h*) 
for  the  Crank-Nicolson  scheme. 

Theorem  4.1. 

1.  If  £  <  p  <  1,  then  Scheme  (4.1)  is  unconditionally  stable. 

2.  IfO  <  p  <  |  and  a  >  0,  then  (4.1)  is  conditionally  stable,  the  stability  condition  is 


r< 


2(1-2  p)(o*  +  6*)’ 


(4.2) 


3.  If 0  <  p  <  |  and  i>0,  then  (4.2)  is  unconditionally  unstable. 
Proof.  The  root  of  the  characteristic  polynomial  for  (4.1)  is 

R  m  1  -  (1  -  -  »(1  -  v)h 

1  +  P07  +  ip* 7 


Therefore 


•  I  “  (1  +  pa7):  +  p:i*7* 


(4.3) 


1' 
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The  condition  |f?|2  <  1  reduces  to 

1  (a8  +  bZ)“l  -  a  <  0.  (4.4) 

If  j  <  p  <  1,  then  (4.4)  automatically  holds;  hence  Scheme  (4.1)  is  unconditionally  stable.  If 
0  <  p  <  then  the  condition  (4.4)  reduces  to  (4.2).  Clearly,  if  a  =  0,  (4.2)  cannot  be  satisfied  for 
r  >  0  and  hence  is  unconditionally  stable. 

I 

Note  that  when  i  *=  0,  (4.2)  reduces  to  the  well-known  stability  condition  for  the  heat  equation 

129]- 

4J.  The  Multi-Dimensional  Case: 

In  this  section,  we  are  going  to  extend  the  results  in  Sec.  4.1  to  the  case  of  multi-dimensions. 
We  suppose  the  aquation  b  of  the  form  (1.2)  and  consider  the  natural  extension  of  the  scheme 
(4.1)  to  multi-dimensions: 

(<S) 

1=1  1=1 


where  we  have  used  the  same  notation  as  in  Sec.  3.3. 

As  in  Sec.  4.1,  here  p  is  a  parameter,  0  <  p  <  1.  Before  we  state  the  stability  results  for  this 
scheme,  we  need  a  few  definitions. 


Deflnttioa  4.1.  Define  an  m-dimensional  index  vector  v  to  be  a  vector  in  IF"  with  components 
having  values  of  either  0  or  1.  Define  Im  to  be  the  set  of  all  m-dimensiona!  index  vectors  except 
the  vector  {0, (),•••, 0}r. 

Theorem  4.2. 


1.  lt\  <  p  <  1,  then  Scheme  (4.S)  is  unconditionally  stable. 

2.  If  0  <  p  <  |,  and  a/  >  0  for  l  *  1, •••  ,m,  then  (4.5)  is  conditionally  stable,  the  necessary  and 
sufficient  condition  is 


k< 


mm 


2(1  -  2p) 


tel  * 


Eth 

tel  ' 


(4.6) 


Proof.  The  root  of  characteristic  polynomial  for  (4.5)  b 


where 


i  -  (i  -  »)  E  Am 


i=i 


i  +  /» E  Am 

tel 


7 1  —  4r;  sin 

k 

fl~  hj’ 


2’ 


0  <  I,  <  2x, 


9 


9 


9 


"•V  «  •  *  •  9\ 
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We  can  easily  obtain 

fl  -  (1  -  p)  £  anil  +  (1  -  m)2  ( E  bni) 

|  pi2  _  l _ <=i  J _ \i=i  / 

-  7  -  rs  i~  xi 

+  m  E  °n*  J  +  /**  (E  *ro  J 

Then,  from  the  stability  condition  jPj  <  1,  we  get 


(4.7) 


tel 


S-,[(S")’*(£"')'h 


(4.8) 


It  is  clear  that  (4.8)  is  true  if  |  <  p  <  1,  and  (4.5)  is  unconditionally  stable  in  that  case.  In  the 
case  0  <  p  <  i,  it  can  be  verified  that  #(7i,”*,7m)  reaches  its  maximum  value  in  the  region 
D  *  {(*»!,•  •  • ,  *y m),  0  <  -n  <  4ri,  /  *  1, • . • , m)  only  at  the  boundary  of  D.  Clearly,  at  (0,  •  •  •  ,0), 
(4.8)  always  holds.  Hence,  we  obtain  (4.6). 


Corollary  4.1.  Suppose  0  <  p  < 

1.  If  ai  =  0  for  some  l,  then  (4.5)  is  unconditionally  unstable. 

2.  If  all  the  bi  have  the  same  sign,  then  the  stability  condition  is 

*-2(1-2?) 

3.  If  bi  =  0,  1  <  /  <  m,  then  the  stability  condition  is 

k<—l _ — 

~  2(1 -2/i) 

teJ  k> 

which  is  a  classical  result  for  the  heat  equation  129}. 


(4.9) 


(4.10) 


Proof.  To  prove  the  corollary,  note  first  that  if  o<  =  0,  then  by  choosing  i>  =*  (0,  •  •  •  ,0, 1,0,  •  ■  •  ,0), 
with  the  *1*  in  the  /-th  position,  in  (4.6),  we  have  k  <  0  and  hence  the  scheme  is  unconditionally 
unstable. 

Second,  if  all  the  b j  have  the  same  sign,  the  minimum  in  (4.6)  must  occur  for  v  =*  (1, 1,  •  ■  • ,  1}T 
from  which  (4.9)  follows. 

Lastly,  if  »  0  ,  1  <  /  <  m,  then  (4.10)  follows  directly  from  (4.9). 

I 

The  result  (4.6)  can  be  viewed  as  the  extension  of  the  stability  result  for  the  multi-dimensional 
heat  equation.  It  is  easy  to  see  that  (4.6)  is  also  the  extension  of  (4.2)  in  Theorem  4.1. 


5.  The  Leap-frog  Schama 


We  now  consider  some  three  level  schemes.  First,  we  study  the  Leap-Frog  Scheme: 


The  truncation  error  is  0(Jt*,  A*). 


(5.1) 
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Theorem  U.  Scheme  ( 5.1 )  is  stable  if  and  only  if  o  *  0  and  r  <  jjjj. 

Proof.  The  characteristic  polynomial  is 

4(z)  ■  s2  +  2 (a  +  *4)7*  -  1.  (5.2) 

We  shall  use  the  Schur  theory  outlined  in  Sec.  2  to  determine  conditions  under  which  4>{z)  is 
Simple  von  Neumann.  We  have 

d*(*)  =  -z2  +  2(a  -  ih)~tz  +  1 
and  4i(t)B4«u 


Since  |^*(0)|  =  |d(0)|,  4>{z)  can  be  simple  von  Neumann  only  if  di(z)  =  0,  i.e.  o  *  0.  Assuming 
this,  for  d(z)  to  be  simple  von  Neumann,  d'(z)  &  2 z  +  2$h  must  be  Schur,  which  leads  to  the 
condition  r  <  jjjj. 

The  stability  result  for  the  Leap-Frog  scheme  can  be  easily  generalized  to  higher  dimensions. 
Consider  the  following  scheme  for  (5.1) 


„•+*  _  u 


2k 


i=i 


(5.3) 


The  corresponding  characteristic  polynomial  is 


*(z)  =  *2  +  2^53Aj7/^  *" 


By  following  the  same  proof  as  outlined  above,  it  is  easy  to  derive  the  following  stability  condition 
for  (5.3) 


10*ni  <  1- 

(=i 

From  (5.4)  it  is  straightforward  to  derive  the  following: 
Theorem  5.2.  Scheme  (5.3)  is  stable  if  and  only  if 

oi  =  0,  l  «  1,2," 

and 

k< 


(5.4) 


,m 


1 


where 


4max{2^,Bw}’ 

*r<0 


t/>0  1  s.<o  1 


6.  The  Du-Fort  Frankcl  Scheme 

This  is  a  well-known  scheme  for  the  heat  equation.  With  regard  to  (5.1),  using  |  (u"+s  +  u") 
instead  of  u"+1,  we  obtain  the  Du-Fort  Scheme: 


—  ||*  A 

- - 1  -  uDF7+'“< 


2k  h V 

where  DF?+,u  s  u"^1  -  ti"  -  u"+*  +  «?+*.  The  truncation  error  is  0  ^Jfc2,  A2,  (g)*)  . 


(6.1) 


•  J 


»  J 


2 


Theorem  6.1.  Scheme  (6.1)  is  unconditionally  stable. 

Proof.  The  corresponding  characteristic  polynomial  is 

4(z)  *  (1  +  2Ar)z*  -  (4 Areozt)z  -  (1  -  2Ar). 

We  thus  have 

4*(z)  m  -(1  -  2(o  -  iljrjz2  -  (4(e  -  H)rco$4)z  +  (1  +  2(a  -  i‘6)r), 

4\(z)  ■  8 art  —  8 ar  cos  #, 

|**(0)|2m(l  +  2ar)2  +  4*V, 
and  |*(0)j2  »  (1  -  2or)2  +  442r*. 

Since  |^*(0)|  >  |$(0)|  and  ^j(z)  is  clearly  simple  von  Neumann,  this  scheme  is  unconditionally 
stable. 

I 

The  extension  of  (6.1)  to  multi-dimensions  is: 

U»+S  _  u«  m 

-L—J-mjrAiDfJ?1*  (6.2) 

Iml 

where  DFfj  is  the  DuFort-FVankel  operator  in  the  zr-direction.  The  characteristic  polynomial  is 


By  following  the  same  proof  as  in  the  one  dimensional  case,  with  ar  replaced  by  k  £  and  4r  by 


one  can  easily  prove: 

i=i  1 


Theorem  6.3.  Scheme  (6.2)  is  unconditionally  stable. 


7.  The  Three  Level  Backward  Difference  Scheme 


This  scheme  [6,  29]  is 


ADj+l«, 


and  its  truncation  error  is  0{k?,h2). 

Theorem  7.1.  Scheme  (7.J)  is  unconditionally  stable. 

Proof.  The  characteristic  polynomial  is 

*(*)  «  (5  +  «*T  +  r2  -  2z  + 

*’(*)  -  \**  ~  2*  +  (|  +  «7  “  , 

•nd  4t (z)  ■  f  Q  +  +  ^7*1  *  -  2(1  +  «if  -  ih)- 


(7.1) 


Page  IS 


It  eao  be  verified  that 

!*’(o)l  >  l*(o)| 

and  therefore  for  stability,  (z)  must  be  simple  von  Neumann,  which  gives  the  condition: 


2(1  +ai-ih)  <  j 

(§  +  «t f)2  +  4*7* 

Let  Te  (§  +  07) 5  4-  42tj,  then  condition  (7.2)  is  equivalent  to 

4  ((1  +  «7)a  +  4*7*)  <  T* 
or  4(r-  5+ «7)  <  r* 
or  (T  -  2)*  +  1  +  4bt  >  0 

which  is  always  satisfied  since  a  >  0.  Therefore,  this  scheme  is  unconditionally  stable. 


The  extension  of  (7.1)  to  multi-dimensions  is 


The  characteristic  polynomial  is 

*(*)  “  «n#  +  •  £  *H(^  **  ~  2z  +  i. 


E  *»!>■ 


(7.2) 


I 


(7.3) 


By  following  the  same  proof  as  in  the  one  dimensional  case,  with  07  replaced  by  on:  and  47  by 

(si 

m 

£  one  can  easily  prove: 

(31 

Theorem  7J.  Scheme  (7.3)  is  unconditionally  stable. 


8.  Schemes  for  the  Real  System 

The  schemes  considered  so  far  are  applied  directly  to  the  equation  (1.1),  which  is  complex¬ 
valued  in  general.  But  note  that  if  we  let  u  =  v  +  iw  (where  v  and  tr  are  real  functions),  then  we 
can  rewrite  (3.2)  into  the  following  real  system: 

vt  -  ovH  -  4w„,  (8.1) 

Wt  *  bv„  +  av„.  (8.2) 

While  any  scheme  for  (1.1)  has  a  direct  analog  for  (8.1)  and  (8.2),  this  new  system  opens 
up  more  possibilities  for  constructing  numerical  schemes  because  the  individual  terms  of  the  right 
hand  side  of  (8.1)  and  (8.2)  can  be  treated  independently  of  one  another  by  different  methods.  It 
is  also  straightforward  to  implement  these  schemes  in  real  arithmetic.  In  the  next  two  sections,  we 
shall  consider  a  few  examples  of  such  schemes. 


8.1.  A  Two-Lsvel  Heal  System  Scheme 

In  this  section,  we  shall  consider  the  following  2-level  scheme  which  is  similar  to,  but  effectively 
different  from,  the  Euler  scheme  for  (1.1). 


Jp  [«  (»"+»  “  2 »"  +  «"-»)  “  *  «+»  ”  2 I®’  +  «"-!)] 

p  [*  («"+?  -  2o"+1 + v"-0 + a  (®T+t  '  2vi  +  «?-i)] 


(8.3) 

(8.4) 


This  scheme  appears  to  be  semi-implicit,  but  actually  it  is  explicit  for  computing  because  we  can 
compute  r"+1  from  (8.3)  explicitly,  then  substituting  v"+1  into  (8.4),  we  can  explicitly  compute 
w*+1.  It  is  easy  to  see  that  the  truncation  error  of  (8.3)  and  (8.4)  is  0(k,h2). 

Theorem  8.1.  Scheme  (8.3)  and  (8.4)  is  stable  if  and  only  if 


1 

2(«+i6ir 


(8.5) 


Proof.  The  amplification  matrix  for  this  scheme  is 


Gm 


(  1  0\~l  (l-ai  h  \ 

\h  l)  \  o  1  -a-r) 

(  i-«t  h  \ 

V-(l-«7)l7  1-07-8*7*/’ 


and  the  characteristic  polynomial  is 


d(z)  *  r*  -  r  (2(1  -  «7)  -  8*7*)  +  (1  -  07)*. 

It  follows  that 

d*(s)  -  (1  -  «7)S*1  -  *  (2(1  -  «7)  -  8*7*)  + 1 
and  *»(*)  -  s  (1  -  (1  -  07)*)  -  (l  -  (1  -  *7)*)  (2(1  -  *7)  -  8*7*)  • 


Since  |4*(0)|  >  |d(0)|,  (z)  has  to  be  simple  von  Neumann  for  the  scheme  to  be  stable. 

Defining 

%  _  2(1  -  87)  ~  6*7* 

l  +  (l-«7)*  ’ 

this  reduces  to  the  condition 

-1  1. 


Since  A  can  be  rewritten  as 


1  +  (1-S7)* 


it  can  be  easily  ieen  that  the  condition  A  £  1  is  always  satisfied.  The  condition  that  -1  «S  A 
reduces  to 


1(7)  m  -(07  -  2)*  +  8*7* , 
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which  b  equivalent  to  the  condition 
and  gives  the  stability  limit 


*r|*|  <  2  -  ay, 

^  1 
*(•  +  W)‘ 


I 


Note  that  this  scheme  spans  a  stencil  identical  to  that  for  the  Euler  scheme  but  the  stability 
condition  is  quite  different.  Unlike  the  Euler  scheme,  thh  scheme  is  conditionally  stable  even  when 
•  a  0.  This  improvement  results  from  the  implicit  treatment  of  the  v-term  in  (8.4)  and  b  a  direct 
consequence  of  separating  the  original  equations  into  real  and  imaginary  parts. 

At  first  sight,  it  may  appear  passible  to  switch  the  role  of  v  and  v  in  (8.1)  and  (8.2)  at 
alternate  time  steps,  similar  in  spirit  to  the  ADI  method,  in  order  to  achieve  a  combined  scheme 
that  is  second  order  in  time.  Unfortunately,  it  can  be  verified  that  thb  b  not  true  and  the  resulting 
alternating  scheme  b  still  first  order  in  time. 

M.  The  Two-Level  Real  System  Scheme  In  Multi-Dimensions 

The  results  in  the  previous  section  can  be  extended  to  higher  dimensions.  Consider  the 
equation  (5.1)  which  can  be  written  into  the  following  real  system 


vi  m  y*!  1gMi  “  5" 'hiwmiMl, 

I*  1  M 

m  m 

w«  *  53  +  XI  •»••!•»* 

tail  tad 


(8.6) 


and  the  following  scheme: 


£e|Z£|t>-£M£|W, 
t*  1  lei 


I«1  1*1 


(8.7) 


The  truncation  error  for  (8.7)  b  again  0(k,  h?). 


m  m 

By  replacing  the  terms  ay  and  5y  in  the  previous  proof  by  J2  at"U  and  £  *171  respectively, 

1*1  tsl 


the  same  proof  goes  through  for  (8.7)  and  we  obtain  the  following  stability  result: 


Theorem  14.  Scheme  (8.7)  is  stable  if  and  only  if 


k< 


84.  A  DuFort-Arankel  Leap-Frog  Scheme  for  the  Real  System 

Although  the  real  system  (8.1)  and  (8.2)  ean  result  in  schemes  that  have  superior  stability 
properties  than  schemes  for  the  complex  equation  (1.1),  deriving  such  schemes  b  not  at  all  au¬ 
tomatic.  In  thb  section,  we  shall  show  that  a  rather  natural  scheme  for  the  real  system  (8.1) 


9 


1...  _ 


» 
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and  (8.2)  actually  has  worse  stability  properties  than  the  corresponding  schemes  for  the  complex 
equation  (1.1). 

We  construct  a  real  3-level  scheme  for  (8.1),  (8.2)  as  follows: 


i  t 


2* 


_  ) 


»— i 


24 


p  (•*« -  >r,-»r,+ >’*-.) -p  «+.  -  *•? + «?-.) .  m 

p  («&■  -  *r + «?.+;) + p  (•?«  -  «r  -  «r  ■ + «•-.)  •  <»•»> 


This  is  a  combination  of  Du- Fort  Frankel  Scheme  and  Leap-frog  Scheme.  The  truncation  error  is 
0(4*,  A*,  ($)*). 

Theorem  U. 

1.  It  •  *  0,  then  (8.8),  (8.9)  (the  scheme  is  equivalent  to  the  complex  Leap-frog  Scheme)  is 
conditionally  stable,  the  stability  condition  is 


(8.10) 


2.  If  I  a  0,  then  (8.8)  (8.9)  (the  scheme  is  equivalent  to  the  Do-Fort  Frankel  Scheme  for  the  heat 
equation )  it  unconditionally  stable. 

3.  It  a  #  0  and  *  ft  0,  (8.8)  (8.9)  is  unconditionally  unstable. 

Proof.  The  Fourier  transform  of  (8.8), (8.9)  can  be  written  in  this  form: 

(1  +  2«r)J5*+1  ■  A/5*  +  (1  -  2er)/5*”1,  (8.11) 


where 


5 


(  «(4r  —  7) 

\  -2h 


2h  \ 

•(4r-*r)y  ’ 


/  is  the  identity  matrix  and  (v,«)  are  the  Fourier  transform  of  v  and  v. 

It  is  easy  to  verify  that  At  is  normal  and  the  eigenvalues  of  At  are: 


p±  ■  a(4r  -  27)  ± 


(8.12) 


ft,-) 


Therefore  At  can  be  written  as  At «  X”1  AX, where  A 
Multiply  from  the  left  of  (8.11)  by  X,  we  obtain 

(1  +  2«7)X5*+i  •  XA/5*  +  (1  -  2a7)XS*”1. 


and  X  b  a  unitary  matrix. 


(8.13) 


Defining 


we  have 


(1  +  2ar)s*+l  -  XAtX~la*  +  (1  -  2«r)s*-‘. 


(8.14) 


Due  to  the  property  of  X  and  the  definition  of  s,  the  stability  of  s  and  u  are  equivalent.  So  we 
need  only  to  discuss  the  stability  of  s.  Since  XMX~l  »  A,  we  obtain  from  (8.14): 


(1  +  2or)s*+I  »  As*  +  (1  -  2or)a*“I 

(8.15) 

(1  +  2«r)y"+I  =  M+y"  +  (1  -  2ar)y*_I 

(8.16) 

(1  +  2«r)*"+1  «  p_r*  +  (1  -  2ar)*"_1. 

(8.17) 

Their  characteristic  polynomials  are 

*±(A)  *  (1  +  2sr)As  -  p±A  -  (1  -  2or). 

(8.18) 

From  this  we  can  easily  derive 

**(A)  «  (-1  +  2ar)A*  -  p,A  +  (1  +  2ar), 

*i(A)  *  1(1  +  2«f)*  -  (1  -  2ar)*]A  -  ((1  -  2 ar)p±  +  (1  -  2ar)p7]  «  0,  (8.19) 

and  d'(A)  *  2(1  4-  2or)A  -  fi±. 

In  the  case  a  «  0,  |d*(0)|  *  |d(0)|,  and  ^t(A)  &  0.  Therefore  stability  requires  ^'(A)  t0  be 
Schur,  which  reduces  to  the  condition  1117  <  1.  The  stability  condition  is  thus  r  < 

In  the  case  a  >  0,  |d*(0)|  >  |^(0)|,  so  ^i(A)  must  be  Simple  von  Neumann.  The  only  roots  of 
*i(A)  are 

Hence 

M-(t -%)**?. 

Clearly,  if  4  m  0,  then  |A±|J  <  1  because  0  <  7  <  4r,  and  (8.8), (8.9)  is  stable.  If  6  /  0, 
|A;t|2  «  1  +  16 42rJ  >  1  when  7  =  4r  and  (8.8), (8.9)  is  unstable. 

I 

Therefore,  as  far  as  stability  properties  are  concerned,  this  scheme  is  similar  to  the  Leap  Frog 
scheme  but  is  inferior  to  the  complex  Du- Fort  FVankel  Scheme,  which  is  unconditionally  stable. 

9.  Summary  and  Concluding  Remarks 

in  this  paper,  we  have  presented  a  rather  exhaustive  collection  of  difference  schemes  for  the 
Schrodinger’s  equation  (1.1).  This  includes  explicit  and  implicit  schemes,  two  and  three  level 
schemes  and  schemes  with  artificial  dissipation.  It  also  includes  schemes  derived  from  the  real 
system  obtained  after  separating  the  solution  into  its  real  and  imaginary  parts. 

While  many  of  the  schemes  are  adaptations  of  well-known  schemes  for  the  wave  equation  and 
the  heat  equation,  the  stability  properties  are  often  quite  different  and  we  have  summarir«d  the 
main  results  in  Table  1. 


Scheme 


Error 


Stability  Condition 


Taylor  Series  for  u*  =  tuM 


0(fc»,fc*) 


p  *  1  :  unstable  (Euler) 
p  -  2  :  unstable  (Lax  Wendroff) 

p«  3:  r  <  ^ 


General  2-level  Scheme 


/  '-.j.  =  p  EAjDff'u 


+(1  -  M)  £ 

<»i 


J*?4  | : 
0(A,AJ) 

OfA’.A*) 


(i)  0  <  p  <  | 

or  *  0  for  some  l:  unstable 

m 

all  1/  have  same  sign:  A  <  gfrig)  •.-■**■ — 

|f$ 

•nh-0:  - 

(ii)  |  <  /i  <  1 :  stable 


Euler  with  dissipative  term  AuMM  (•  ■  0) 

+  (o|6|  +  ifib)k(Dj)2* 

Scheme  with  least  restrictive  stability 
*  ibDjU  +  (-£  +  f )  k(Df)*u 


0(k,k 2) 


0[k,h2) 


(a  <  P  >  0  or  a  <  — |t  l  <  0) 

t+^/-a(2o*+29*+a) 
f  ^  ~  4(a»-4>)W - 


r$5W 


Euler  with  dissipative  term  A*urrr,  (a/  =  0) 


a  <  0 
and 


+  (° + »vj)a?  (f^,)’*) 


<«i 


Scheme  with  least  restrictive  stability 

('*>>  +  !(-'+  0  '>!(^)!“) 

tel  ' 


0(A,A») 


0(A,A*) 


A  <  min 


v  1-1  • 


*s-J— 

tail  1 


Table  I:  Summary  of  Schemes  for  u t 


1  22  AiumiBl,  Ai  =  o/  +  »'6|,oj  >  0,  #■#  + tw. 

tel 


Next,  we  would  like  to  comment  briefly  on  the  effect  of  lower  order  terms  (such  as  u,  and  w) 
on  the  stability  results,  which  arise  in  many  applications.  If  we  adopt  the  slightly  weaker  definition 
of  stability  (2.1),  under  which  the  necessary  and  sufficient  stability  conditions  given  in  Table  1  are 
still  sufficient,  then  it  can  be  shown  [29]  that,  when  lower  order  terms  are  present,  these  conditions 
remain  sufficient  to  insure  convergence  of  the  numerical  solution  in  the  limit  as  A  and  A  tend  to 
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Scheme 

Error 

Stability  Condition 

Leap-Frog 

Of*2,  A*) 

all  a/  =  0  and 

b  —  41 UX(1,,£JV}’ 

where  #Jv  =  £ 

*,x>  1  b,< o  *< 

Du-fort  Frankel 

tel 

0  ({)*) 

stable 

Backward  Difference 

of*2,  a2) 

stable 

2-level  Real  System  Scheme 

HI  ffll 

J  '  i  m  f1  fl.n»  «_  f*  A.m 

0(4,  A*) 

k "  1 

l  *  a‘L/iJtv 

„«+i  *n  m 

-  £  hoZ" + £ 

tel  tel 

Du-fort  Frankel-Leap  Frog 

Real  System  Scheme 

i 

'  m  aDF7v  “  bDlw 

•41  »•! 

=  6DJ+Io  +  aDFJv 

•  *0:  r<^, 

6*0:  stable 

«  /  0  and  6  j*  0 :  unstable 

Tabls  I:  Continued 


sero  (assuming  the  scheme  is  consistent).  However,  for  fixed  k  and  h,  the  numerical  solution  may 
grow  with  each  time  step  while  the  exact  solution  decays  with  time.  We  plan  to  further  investigate 
the  practical  stability  of  these  schemes  in  the  presence  of  lower  order  terms. 

Stability  results  are  important  for  at  least  two  reasons:  to  ensure  convergence  of  the  numerical 
solution  to  the  exact  solution  as  k  and  h  tend  to  sero,  and  to  ensure  that  the  numerical  scheme  is 
insensitive  to  round-off  errors  for  fixed  k  and  h.  The  choice  of  a  numerical  scheme  for  a  particular 
application  depends  on  more  than  just  the  stability  property.  In  practice,  accuracy  and  efficiency 
are  often  the  controlling  factors.  The  optimal  choice  can  often  only  be  detemined  empirically  with 
some  experimentation. 
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While  we  have  limited  ourselves  to  second  order  centered  finite  difference  schemes  for  the  spa¬ 
tial  discretisation,  and  have  not  considered  the  use  of  higher  order  and  finite  element  and  spectral 
methods,  the  framework  used  here  can  be  used  to  derive  thses  schemes  and  their  corresponding 
stability  results  as  well. 

10.  Appendix:  Proof  of  Theorem  3J 
Proof.  The  amplification  factor  is  given  by 

«  ■  1  -  ^  k/fi  [stw  -  (o  +  iflW] 


where 


k 

rt~  hf’ 

.  .  jfl 
m  *  <  sin* 


Stability  requires 


1*1*  *  [  1  +  £  ai/W]  +  TjT(k/r/,»  -  fibiurft)  <  1. 

I  lal  i  ilatl 

w  /  m  \  *  r  m  1* 

2  J^okjnqf  +  a7  I  i«|f  ]  +  £(k|f|i|l  -  0kfl»l<)|  <  0. 
Isl  Vial  /  Ll«l  J 


P  " 


If  we  define 


then  this  condition  can  be  written  as 


1 53  W  +  a*k(^2fttif)  +  *  ~  £0, 

1*1  Vial  /  Vial  lal  / 


2a  £  /inf 

lal 


/  m  \*  /  »  m  \  2 ' 

*(£*)  ♦('£**  *2**) 


P®  (?!.•••  >»Jm) 

£“{P*»  0  <  in  £  4, 1  «!,•••, m) 

{«  m  \ 

p;  p€Z?  and  £  fm  <  fi £ /i«|f  > 
la!  lal  J 

0**(p;  p€D  and  ■ 

l  lal  In  ] 


and  define 
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Clearly,  D  *  Di  U  Di- 

We  separate  the  two  eases:  0  >  $  and  0  < 

Cast  1.  0>\ 

1.  In  Di 

Since  _  _ 


we  have 


52//H  -  fiUto*  “d  0  -  ^  -  4- 
o  <  -^2  fm  ^  p  52  fM*  ~  \  Ylf,r,t 


and  therefore  „  _  ,  , 

r  ^  _ 2q  l.  fwf _ 

2a 

a 

“-8  [a’ +  (/*-!)*]£/»' 

On  the  other  band,  in  the  ease  0  >  \,  (4,--*,4)  €  A.  and 


G  >  - 


<7(4,--,4)»-; 


RTFrijtf^- 


Hence 


2.  In  Ih 
Since 


we  have 


mf  8[a’  +  (£-J)*]E//' 

E*»  *  (DO*  (£**)'• 

o<Eto-/>E/«?  -  (E^)'  (Eto)’  -  ^Eto*. 


and  therefore 


Letting 


we  have 


G  >  - 


«.*(£  //I?)  +  [(£/)*-  />E  to?)’]  ’ 


G  >  S(w)  s  - 


.(Eto?)1. 

_ 2a _ 

o*wa  +  [(£/iM  “  Pw] 


(10.1) 


fl 
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It  is  easy  to  see 


inf  •  *  0, 

Ol 


and  if  we  define 


w*  «  sup  », 

D , 


then  dearly, 


•*  £  w*e>(4,---,4)«4(£//)J. 

Because  of  the  form  of  S(w),  it  achieves  its  minimum  at  the  end  points  of  the  range  of  v  and 

G  >  min(5(0),  S(w*)) 


Hence 


igf  •  -.ij*)  -  min  S(«’)^  • 


(10.2) 


From  (10.1),(10.2),  we  have 


* s  ""  ('JFTi^TFTO'  '&•  s<”'))  • 


S(m*)  >  S(w)  s -■ 


we  have  the  following  stability  condition  when  /f  >  4, 


* s  mi°  ('«[?■+(#- D*iE/r  ~E7i)  • 


(10.3) 


Casa  JL0<1 

Clearly,  Dt  is  empty  and  D  •  D%.  It  b  easy  to  see  that 


G  >  min  (s(0),S  (<  (£»*)) 


The  general  stahiUty  condition  b  therefore  (10.3).  Clearly,  we  must  have  o  <  0.  To  choose  o  and 
$  such  that  the  stability  condition  b  the  least  restrictive,  we  should  take  $  *  \  so  that  we  have 


tSni,,('5T^’'E7i)- 


To  maximise  the  right  hand  side,  we  should  take 


—  "A 


and  the  stability  condition  becomes 
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